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ABSTRACT 



For a relativistic particle under a constant force and a linear velocity dissi- 
pation force, a constant of motion is found. Problems are shown for getting the 
Hamiltonian of this system. Thus, the quantization of this system is carried 
out through the constant of motion and using the quantization on the velocity 
variable. The dissipative relativistic quantum bouncer is outlined within this 
quantization approach. 
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1. INTRODUCTION 



It is known that dynamical systems with dissipation (dissipative systems) 
represent some difíiculties for their right formulation in terms of Lagrangian 
and Hamiltonian formalisms (López, 1996) . Normally this dissipation is in- 
cluded in the dynamical equations in a phenomenological way and through a 
force which depends on the velocity of the particle. If the Hamiltonian of a dis- 
sipative system is found, one proceeds to try to quantize the system. This has 
been able to do for several nonrelativistic systems (Glauber and Manko, 1984; 
Okubo, 1981; López and Gonzàlez, 2005; Mijatovic et al, 1985), but little is 
know about relativistic dissipative systems. In this work, a constant of motion 
for a relativistic particle under a constant force and a linear dissipative force 
is obtained. The constant of motion at first order in the dissipation parameter 
and the nonrelativistic limit of the constant of motion are analyzed, and the 
problems for getting the Hamiltonian of the system are outlined. Finally, the 
quantization for the relativistic dissipative system, at first order on the dissi- 
pation parameter, is carried out through the quantization of the velocity and 
the constant of motion. 



2. RELATIVISTIC CONSTANT OF MOTION 

The one-dimensional motion of a particle of mass m at rest under a constant 
force, /, and a linear dissipation force, — av, is governed by the equation 




mv 



1 — v 2 /c 2 



if + av), (1) 



where v is the velocity of the particle, a is the dissipation parameter, and c is 
the speed of light. This equation can be written as the following autonomous 
dynamical system 

dx , n s 

- = „ (2a) 

+ (26) 

where (3 is the constant defined as j3 — a/ f . A constant of motion for this sys- 
tem is a function K = K(x,v) which satisfies the following partial differential 
equation of first order (López, 1999) 

v 9 -f - kl + Pv)(l - v'/cY 2 ^- = • (3) 
ax m av 

The general solution of this equation is given by 

K p (x,v) = G(A(v)+jx) , (4) 
where G is an arbitrary function, and A(v) has been defined as 

A{v) = m J (l+Pv)(l-v*/c*)s/* ■ (5) 
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The result of the integration of (5) is 



A(v) 



[1 — (3v)mc 

(j)(v 



+ 



mj3c 3 



. dc + v/c 
arcsm -, . V 
1 + pv 



(1 - (3 2 c 2 ) 3 / 2 

mc 2 (l-2v/c-2v 2 /c 2 ) 

1 - v 2 /c 2 ' 3{l + v/c)\Jl-v 2 /c 2 



mcv 



+ 



(1 — (3v)mc 2 



+ 



mj3c 3 



2„2 



log 



if (3c < 1 
i//3c=l 
if (3c>\ 



(6) 



where the functions <p(v) and ip{y) have been defined as 

(f>(v) = yjl - V 2 /C 2 (1 - /3 2 C 2 ) 

and 

ij)(v) = 2 (^c 2 /3 2 + (3v + yjl - v 2 /c 2 ^Jc 2 (3 2 - 1 



(7a) 



(76) 



The function G, appearing on (4), can be determinate through the criterion 
(López, 1996) of having the usual relativistic energy expression for (5 equal to 
zero (non dissipative case), 



mc 



lim Kg(x, v) = 



+ fx, 



(8) 



which brings about the result G = I (the identity function) 1 . Therefore, the 
constant of motion for the system (2) can be chosen as 



K p {x,v) = A{v)+fx . 



(9) 



This constant of motion brings about the damping effect on the trajectòries in 
the phase space (x,v). Of course, due to multivalue functions of (6), the value 
of the constant of motion changes for the trajectòries going from the x > 
side to the x < side of this space, in order to get the spiral falling down to 
the origin behavior of the trajectòries. So, one may say that (9) represents an 
"almost everywhere" constant of motion of the system (2), in the sense that 
the set of points where these changes do occur has zero measurement (Hewitt 
and Stromberg, 1965), one one may call it "local constant of motion." 

For weak dissipation, one can also make a Taylor expansion on (5) of the 
term (1 + f3v ) _1 to get the constant of motion of the following form 



K p {x,v) = 



mc 



1 — v 2 /c 2 



■m(3 



vc 



yjl - V 2 /c 2 



c arcsm - 

c 



+ fx + $(v) , (10a) 



1 one could add to (8) the term —mc 2 to have the usual energy expression for the nonrelativistic 
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where the function $ is given by the expression 



n=2 



c 2 v n 



nc 



n-l)Jl-v 2 /c 2 



n 



v n 1 dv 



(l-v 2 /c 2 ) 3 / 2 



(106) 



Thus, at first order on the dissipation parameter (3, one has 



Kp(x, v) = 'jmc 2 — m(3c 3 
where 7 has the usual expression, 
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71» . v 
arcsm - 

c c 



+ fx 



1 



y/l ~ V 2 /C 2 



(lla) 



(116) 



Note that the nonrelativistic limit {v/c <C 1) must no be taken from the 
case C/9 < 1 or the case cf3 = 1 but rather from the case cf3 > 1. In this 
way, subtracting the term of energy at rest, mc 2 , on the case c/3 > 1, the 
nonrelativistic constant of motion is given by 



v í \ m 



j3v - log (1 + pv) 



+ fx. 



(12) 



Now, using the known expression (López and Hernàndez, 1989; Kobussen,1979; 
Leubner, 1981) to obtain the Lagrangian from the constant of motion, 



K(x, v) dv 



(13) 



one can get the Lagrangian and generalized linear momentum {p = dL/dv) 
given by 

Lp(x,v) — B(v) — fx (14a) 

and 

p(v) = C(v) , (146) 

where the functions B(v) and C(v) are given in the appendix. in particular, 
the Lagrangian and the generalized linear momentum for a relativistic particle 
with dissipation at first order in the dissipation parameter f3 (expression lla) 
are 



L^x.v) 



mc 



7 



— m/3c 3 



v , 7 — 1 arcsm ^ 
— log 1 ç - 

2c 7 + I 7 



v 1 
+ - log 



and 



p(v) = •ymv + f3mc 2 



5 v 7 . v 3 v 

1 arcsm log - 

2 c c 2 6 c 



— fx (15a) 



(156) 



Similarly, for the nonrelativistic case with dissipation (expression (12)), one 
has 



Tíí 

L 2 (x,v) = -(f3v-l)\og(l + Pv) , (16a) 
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and 

p(v) = — |log (1 + f3v) 



(166) 



l + /3v. 

As one can see from (14b), (A 2 ), (15b) and (16b), it is not possible to have the 
inverse relation v = v{p). Therefore, their Hamiltonians are expressed only 
in an implicit way through the constants of motion (9), (12) and (11). Thus, 
the quantization of the system (1) can not be carried out with the Standard 
Shròdinger equation (Messiah, 1958), 



or Heisenberg equation, 



ih— = H(x,p)V, (17a) 



4 = < 176 > 



where £ is any time depending observable. The Feynman path quantization 
(Feynman and Hibbs, 1965) is, in principle, possible to use here since one 
has gotten the Lagrangians (14a), (15a) and (16a). However, the analyti- 
cal functions appearing in these expression represent a real challenge for the 
quantization with the path integration method. One must also note from rela- 
tions (14a) and (15a) that these can not be expressed in a covar iant way since 
Lorentz transformations do not leave invariant this dissipation system. 



3. QUANTIZATION OF THE CONSTANT OF MOTION 

We are interested here in the quantization of the system (1) at first order 
in the dissipation parameter f3, characterized by the constant of motion (lla). 
As it was mentioned above, the quantization using the Hamiltonian or the 
Lagrangian approaches does not look plausible due to the implicit form of the 
Hamiltonian and the complicated expression for the Lagrangian. However, one 
can try to use the idea of quantizing the velocity (López, 2000) through the 
obvious expression 

. h d 

mdx 

In this way, one can use directly the constant of motion of our system to make 
its quantization through the equivalent Schròdinger equation 

ih— = K(x,v)V, (19) 

where K is the Hermitian linear operator associated to the constant of motion 
K and which has units of energy. 
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Thus, let us consider the constant of motion (lla). From (19) and (lla), 
the equation obtained is given by 



rac 



V 1 - ^ 2 A 



- ra(3c 3 



v 

— arcsm - 



+ fx}^ 



(20) 



Equation (20) represents an stationary problem, and the proposition 

E V 



ty(x, t) = ïp(x)exp 
transforms (20) to an eigenvalue problem, 



rac 



1 — v 2 /c 2 



— raf3c 3 



V 1 " d '/ c 



— arcsm 



+ fx } ip = Etp 



(21) 



(22) 



One can see from this expression that it is better to look for its solution in the 
velocity representation, which is given by applying the Fourier transformation 
to the function ip, 



(f)( v ) = F[tp(x) 



— le> 
J2ïi Jr 



imvx/h 



ip(x) dx , 



(23) 



where the variable v represents the velocity of the particle. Applying this 
Fourier transformation to (22), one gets 



rac 



1 — v 2 /c 2 



ra/3c 3 



= — arcsm - 

s. c 



+ i 



h£d_ 

ra dv 



= E(p . (24) 



This equation has the following solution 



<Pe{v) 



2 S 

,- m c 



arcsin Zftcu 1 



f 2 v Ev 
— — pv arcsin o 

c z c rac 



(25) 



where v represent some maximum velocity of the particle (v < c), and 
[—l'o, v ] is the intervale of velocities where the normalization of the function 
(25) has been carried out, 



/Vo 
\M 
-Vo 



v)\ 2 dv = 1 . 



One also has that 

{4>e\4>e' 



U* E {v)<t> EI {v) dv=^-5(E-E>) 
J ra 



(26) 



(27) 



The spectrum of energies of the particle is continuous because of the form of 
the potential, fx, and the general solution of (20) can be written, using (21) 
and (25), as 

$(i;,í) = I A(E)Mv)e- lEt/h dE , (28) 
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where the coeficient A(E) is determinate by the initial condition <&(v, 0) in 
the following way 

A(E) = ^j ( j>* E (vMv,0)dv. (29) 

Now, in the case the potential be of the form 

!fx if x > 
(30) 
oo if x < 

which would correspond to the one-dimensional dissipative relativistic bouncer 
problem, one requires that ip(0) = and ip(x) = if x < for the solution 
of (22). This condition brings about the discrete spectrum of the system 
and can be written in the velocity representation, using the inverse Fourier 
transformation, in the following way 

= ^(0) = F-'iMv)} = 4= / Mv) dv . (31) 

x=n \ l'K J 



Given the set of eigenvalues, {E n }, the eigenfunctions are given by 



■ v on l·i 1,2 a ■ v EnV 
arcsin 2pcy 1 — pv arcsin , 

c V c 2 c mc 



(32) 



and the general solution can be written as 

<S>(v,t)=J2 A nMv)e~ iEnt/h , (33) 

n 

where, using the orthogonality (0 n |0 n ') = <5 ni „/, the coefficients A' n s are deter- 
minate by the initial condition $(v, 0) and through the following expression 

An= i" <j>* n (v)$(v,0) dv . (34) 

J — v 

The wave function in the x representation is gotten through the inverse Fourier 
transformation, 

i, n {x) = 4= / e~ mvxln K{v) dv . (35) 



. 2tt 

Note that the above quantization process can be done in general for the con- 
stant of motion (9) and any case defined by (6). 
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CONCLUSION 



For a relativistic particle under a constant force and a dissipative force pro- 
portional to the velocity, a local constant of motion has been found , and one 
has outlined the problem to get its Lagrangian and Hamiltonian in general. To 
quantize this system , Schròdinger quantization approach has been used with 
this local constant of motion at first order on the dissipation parameter and 
using the velocity representation of the wave function. Using this quantiza- 
tion approach, we have outlined the dissipative relativistic quantum bouncer 
problem. 
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APPENDIX 



The function B(y) is given by 



mc 



T3 



1 - 



2 i+A/ 1 -^ 

i-4 + /^iog V / 

c 2 ^ f v j c 



+ (i_ /3V) 3 / 2 ^ 



'2{\ + j\-v*/<? )\ mc 2 f 2 (v) 
B(v) = { -mcvlog ( + 



if /3c = 1 



mc 



-Jl-S + ^log- 



i -/f 2 ? 

where /i,/2 and /2 are defined as 



v/c 



m6c 3 f 3 (v) . £ o i 



(1 _ ^2)3/2 



fi(v)=v j 



. f (3c + v/c \ 

arcsm — av 

l 1 + (3v I 



h{v) 



v v 2 3v I ^, (2(1 + Jl - v 2 /c 2 ) 
-1 h — H \ll log I — 

v/c 



+ - + 

c c z c 



and 



Jl-v 2 /c 2 



(Ai) 



(«2) 



log 



l + Pv 



dv 



fs(v) = v 

where the function ip is given by (7b). The function C(v) is given by 



(«3) 



C(v) = 



mc 



v/c 



—mc 



log 



2(l + Jl-v 2 /c 2 )\ 1 



v/c 



Jl - v 2 /c\ 



+ (1 _ 02^/2 * PC<L 



+ m£9M i{(3c=1 



mc 



1-fF? 



if /3c>l 



(1-/3V) 3 ' 2 



(40 
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where g±, gi and g% are defined as 



9i(v) 
92(v) = 



. fPc + v/c 
f ( \ arcsm . // 



+ 

V V 
921 (v) + 922(v) 



and 



(1 - v 2 /c 2 f/ 2 {\ + y/l - v 2 /c 2 ) ' 

= + - log . 

v v \1 + pv J 



The functions #21 and #22 have been defined as 



v 3v 2 v 3 



sa(B) = ( -4- + )(1 + Wi_ ) 



and 



#22^) = 3(1 - -=■) 



V 2 , / V 2 



log 



'2(1 + ^jl-v 2 /c 2 ) 

v/c 



(Pl) 
(*l) 
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